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I. INTRODUCTION 



In recent years, it was realized that chaos can have useful applications. Among others we mention the diffusive 
ionization of a Hydrogen atom by a microwave field (l]-|3| , the dissociation of polyatomic molecules by laser radiation 
due to the transition to chaos , the stabilization of chaotic laser radiation |j) by the controlling chaos method , 
secret communications Q using synchronization of chaotic systems |J , and chaos is even promisible in the generation 
of music variations and other sequences of context-dependent symbols [l0|| . All applications mentioned dealt with 
classical chaos. 

Recently the main activity in chaos has been shifting to the studies of quantum chaos, particularities in the 
behavior of quantum systems that are chaotic in the classical limit fll|-|l6f|. As a rule, quantum effects distort or 
suppress "useful" manifestations of classical chaos. For example, in the most studied system, a Hydrogen atom in 
a microwave field, quantum effects suppress diffusive ionization by the mechanism of quantum localizatio n |p"7|| (the 
analog of Anderson localization in a solid state jll|) or by the appearance of "scars" of a wave function |L9| in the 
vicinity of classical unstable periodic orbits of the corresponding classical model of the atom |2(J . The only probable 
useful application known to us (and in particular in the area of quantum optics) for quantum chaos is the suggestion 
pl| to utilize the effect of wave function localization in a one-dimensional quantum chaotic system driven by a periodic 
external field, for the generation of electromagnetic field in the Fock state. 

In this paper we discuss the possibility of squeezed states of the light generation p^-|25|] at the transition to quantum 
chaos. We show that the maximal possible degree of squeezing, achievable during some time interval of measurement, 
for chaotic dynamics is much greater than for regular dynamics during the same time interval. We find a direct 
correlation between the degree of squeezing and the value of the local Lyapunov instability rate in the corresponding 
classical chaotic system. 

Our consideration |2(| is based on three simple but general ideas: 

(i) A free electromagnetic field is in a coherent state that is a Gaussian wave packet. A wave packet spreads when it 
propagates through a nonlinear medium. However, there still exists a time interval of well-defined quantum-classical 
correspondence during which the wave packet follows a path in phase space closed to the path defined by (semi)classical 
equations of motion. During this time interval, localized wave packet demonstrates squeezing in one direction of phase 
space and stretching in another direction. As a result, one of the field quadrature components may be observed in 
the squeezed state at some time moments. Due to the presence of the local strong (exponential) instability inherent 
in the underlying classical chaotic dynamics [j3l| , the stretching and squeezing of the wave packet for quantum chaos 
is much stronger than for the case of the regular and stable dynamics, when the distance between two initially closed 
trajectories in phase space increases in time only in a power wise way |plf . 

(ii) The time interval of the quantum-classical correspondence depends on the number of quanta involved in the 
nonlinear dynamics: it has a power dependence on the number of quanta for regular dynamics, and this dependence 
is logarithmic for chaotic dynamics p2 11 . Tq] . In spite of the time interval of the localized wave packet motion, 
when enhanced squeezing is expected, at chaotic dynamics may be very short, this time scale seems quite observable 
in modern experiments on light squeezing [p2|-p5[, where average number of photons N involved in the nonlinear 
interaction is as great as 10 3 — 10 12 . Moreover, in the conditions of the experiments, the time of the localized wave 
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packet motion and the well-defined quantum-classical correspondence may be of the same order as the time scale 
when dissipation and other factors restricting squeezing could be neglected. 

(iii) Along with stretching and squeezing, another inherent characteristic of chaos in bounded phase space is folding. 
For strong chaos there are strong and multiple foldings of the wave packet; therefore, the time intervals of squeezing 
become very short and squeezing even becomes unstable 28 1. However, at the same time that strong folding appears, 
the wave packet becomes delocalized and well-defined quantum-classical correspondence, which we understand here as 
the motion of a wave packet center along a classical trajectory, is broken |33| ]. Therefore, if one studies the dynamics of 
the wave packet only on a time scale of well-defined quantum-classical correspondence, the squeezing could be strong 
and still rather stable for not very strong chaos or even for regular trajectories lying in phase space near the chaotic 
region. 

It is already known that light squeezing can be increased near the bifurcation points between different dynamical 
regimes j$4 25 1 . Among with well-studied parametric media |54| , such an increase of squeezing was predicted also for 
the interaction of the field with 2-level atoms inside a high-Q cavity |55|,[56| . The explanation of enhanced squeezing 
near the bifurcation point has been used |37| and it is very similar to the arguments presented above that utilize 
quantum chaos for large squeezing: At the transition (bifurcation) between different dynamic regimes, there should 
exists some diverging variable, and then the conjugate variable should be strongly squeezed. In spite of the similarity, 
our suggestion has two significant differences. First, all papers devoted to the study of the enhanced squeezing in the 
vicinity of instability threshold p5|]34| |3C| dealt with the integrable or near-integrable systems with regular dynamics. 
Instead, we suggest the use of the transition from regular to chaotic dynamics. Second, from the general viewpoint, 
the transition to chaos and chaos itself are frequent phenomena and, as a rule, they take place for a rather wide region 
of the control parameters, in contrast to a commonly narrow region of control parameter characterizing a bifurcation 
between different regimes of regular motion. 

In this paper we consider the squeezed light generation by nonlinear nonintegrable optical systems obeying the 
transition from regular to chaotic dynamics in the classical limit. As for other problems of quantum chaos [pl]-^6[, 
we deal with the semiclassical limit when a great number of quantum levels N 3> 1 are involved into the dynamics. 
Our consideration is valid for any nondissipative quantum system with one and half degrees of freedom, but we 
demonstrate our main results on the enhanced squeezing for some particular model of the nonintegrable optical 
system: the nonlinear oscillator periodically forced by the classical field. To investigate the nonlinear dynamics of 
the systems and the dynamics of squeezing in the semiclassical limit, it is natural to use the cumulant expansion 
technique psf as a variant of the general 1/iV-expansion method |39[| . In this paper, we use the 1/iV-expansion 
method suggested in |36[ ], which is well-adopted for the problems of quantum optics. We show that as long as the 
wave packet is localized, its dynamics may be well described by the behavior of the mean values accounting for 
quantum corrections and the lowest-order quantum cumulants. We demonstrate that the equations of motion for the 
second order cumulants are essentially the same as those are used in the definition of the Lyapunov exponent for the 
corresponding classical system. This fact gives us the opportunity to find direct correlation between the degree of light 
squeezing and the degree of local instability in the corresponding classical system. We also compare our approach 
with the generalized Gaussian approximation |2^ , |38|| , which is used to find approximate solutions of many problems 
of quantum optics p0|-^2]| , as well as with a similar variant of the cumulant expansion f43j] , and find good agreement 
with minor differences. 

Our presentation is organized as follows. In the next section, using l/W-expansion method and starting from 
the Hciscnbcrg equations of motion, we derive the self-consistent set of equations for the mean values and second- 
order cumulants describing the dynamics of quantum fluctuations in the semiclassical limit for an arbitrary quantum 
system with one and half degrees of freedom. We show that these equations coincide with the equations used for the 
calculation of the classical Lyapunov exponent. In section 3, we discuss the conditions of validity of our basic equations 
for regular and chaotic dynamics and compare our approach with the results of other semiclassical techniques used 
in both quantum optics and quantum chaos studies. We consider the dynamics of light squeezing at the transition 
to quantum chaos, and compare it with the dynamics of squeezing for regular and stable motion in section 4. We 
illustrate our results on enhanced squeezing at the transition to chaos in the model of a kicked quantum oscillator in 
section 5. Our conclusions arc summarized in the final section 6. 



II. SEMICLASSICAL DYNAMICS OF QUANTUM FLUCTUATIONS: GENERAL FORMALISM 

We begin with a single-mode quantum system described by the Hamiltonian H (b, , t) including an explicit depen- 
dence on time, where b and b^ are Bose operators ([&, 6^] = 1) of the mode. We use normal ordering of operators. Let 
our system include some large parameter N 3> 1. The parameter N may represent the number of quanta (photons) 
pumped to the system pfl] or the number of degrees of freedom of a quantum system H . Because we are interested 
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in semiclassical limit N — > oo, it is useful to introduce new operators for the annihilation and creation of photons 

a = b/N 1 ' 2 , ^ = (1) 

with the commutation relation 

[a >a t] = l/JV. (2) 

In the classical limit (N — > oo), one has two commuting c-numbers. The natural quantum states for the consideration 
of the quantum system in semiclassical limit are coherent states [fll| , [39[ . Thus we suppose that our quantum system 
is initially in the coherent state \a) = exp(Naa^ — Na*a)\0) corresponding to the mean number of quanta ~ N. 
Following the general scheme of the 1 /W-expansion method ]39|,[56| , we rewrite the Hamiltonian H in the form 

H = NH N (a,a\t), (3) 

where the Hamiltonian Hm generates the correct classical equations of motion in the classical limit N — ► oo. In the 
classical limit, the operators a, aJ, following (|^), may be considered as c-numbers of order of unity and all quantum 
corrections being of order of or less than 1 /N could be treated by perturbation theory. 

In order to illustrate the representation of Hamiltonian in the form (||) , consider a particular model of the noninte- 
grable quantum system: the nonlinear oscillator periodically driven by the classical field. In the interaction picture, 
the Hamiltonian has the form (% = 1) 

H = Atfb + |& t2 6 2 + sN l ' 2 {b + tf)F(t), (4) 

where b and describe a single mode of the quantum field and k is proportional to the third-order nonlinear 
susceptibility of a nonlinear medium. The last term in (||) corresponds to a coupling of the oscillator with an external 
classical periodically modulated field containing a large number of photons N ^> 1 (e is a coupling constant, F(t) is a 
periodic function of time, and A is a detuning of the mode frequency from the carrier frequency of the external field). 
Using new operators a, a? given in ([!]), the Hamiltonian (Q) may be represented in the form (j^) with 

H N = Aa^a+^a^a 2 + e(a + a f )F(t), g = kN. (5) 

It may be shown that in the classical limit J\f»l, when we have classical variables a, a* instead of operators a, a' 
with | a |~ 1, the dependence g ~ N correctly performs the time scale of the energy oscillation for Kerr nonlinearity 
Q| . We will return to the study of the nonlinear oscillator (^) in section 5. 

We now turn to the general case (^) and derive the equation of motion for mean values and first-order cumulants. 
From the Heisenberg equations for a, a 2 and their Hermitian conjugated equations, we have the following equations 
of motion for the averages over coherent states 

l Tt = (n (6a) 

i ±((6a) 2 )=2(V6a) + (W), (6b) 

i— (Sa*Sa) = -(V*5a) + (5a* V), (6c) 
dt 

where we have introduced 

V(a,a*) = ^-, W(a,a*) = ±^, (7) 
oa> N oa' 

and z = (a), 5a = a — z, ((5a) 2 ) = (a 2 ) — (a) 2 , (5a*5a) = (a^a) — (a^)(a). In the derivation of equations (Q) we used 
the equalities Q 

[a,M(o,at)] = [M(o,at),ot] = - — , (8) 
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which are valid for an arbitrary function M of operators a, (g). 

The set of equations (^) is not closed and actually is equivalent to the infinite dynamical hierarchy system for 
moments and cumulants. To truncate it we make the substitution a — > z + 8a, where at least initially the mean z ~ 1 
and the quantum correction \Sa(t = 0)| ~ N~ x / 2 <C 1. We expand functions V(a, a*) and W(a, a*) around the mean 
value z = (a), 



Sa* + ■■■, (9) 



where the subscript z means that the values of V, W, and their derivatives are calculated at mean values z and z*. 
Substituting expansions (^) into Eqs. (^|) and taking into account the equality 



/dV\ _ I d 2 H N \ /dV_ 



\ da / \ dada* / \da 



(10) 



resulting from Eqs. (Q), we have in the first order of 1/N the self-consistent set of equations for mean values and first 
order cumulants 

ij t z = (V) z + g[z, z\ ((Sa) 2 }, (5a*6a)}, (11a) 

4^ - 2 ^ + 2 {w). {5a * 5a) + {w) - (llb) 

,|<^ fa > = «*«)'> + (f£) «*•)">. die) 



The small quantum correction q ^ 1/N involved in the equation ( 11a ) has the form of the second differential of V as 
follows 



The initial conditions for the system ( pi] ) are ((Sa) 2 )(t = 0) = (6a*Sa)(t = 0) = and some arbitrary z(0) is given, 
which is of the order of unity. 

We now turn to the classical equations of motion. The classical limit may be obtained from (11a) by neglecting the 
quantum correction q of order 1/N and then the classical Hamiltonian equations are 

Linearization of classical equations ( |l3| ) near z by means of the substitution z — > 2 + Aa (| Aa| <C |z|) gives 

d <9F 97 d <9V* 3F* 

i— Aa = — Aa+ — Aa*, i— Aa* = — Aa-- — Aa*, (14) 

at az az* at az az* 

where all derivatives are taken on the classical trajectory found from the Hamiltonian equations (|l3|). Using the 
classical analog of Eq. jl0|), 

dV dV* 



dz dz* ' 

we have from ( |l4j ) the following equations of motion for quadratic variables (Aa) 2 and |Aa| 2 

,|(Aa)' = 2f (A«f + 2 |r|A„P, (15.) 

4|Aa|* = -^l(Aaf + §(Aa-f, (15b) 
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with the initial conditions 



(Aa) 2 (0) = [Aa(0)] 2 , |Aa| 2 (0) = [Aa(0)] x [Aa*(0)], 



(16) 



where Aa(0) and Aa*(0) are the initial values of small deviations from the classical trajectory. We define the distance 
D c i in classical phase space between two initially closed trajectories as 

1/2 



D c i(t) = \Aa(t)\ = (ReAa) z + (ImAa)' 



(17) 



where the time-dependent quantity |Aa(t)| is determined from linearized classical equations of motion (|l4| ) or from 
equivalent equations of motion for the square of linear deviations (|l^) . The value D c i characterizes the degree of local 
instability in the classical system: For chaotic motion D c i(t) increases on average exponentially with time and for 
regular motion D c i has only a power wise time-dependence j3l| . Using D c i(t), one can define the largest Lyapunov 
exponent as 

In Del (*) 



A = lim 



t 



(18) 



If A > 0, the motion is chaotic and A = for regular motion |31[ . 

We turn again to the quantum dynamics. It is easy to see that the equations of motion ( lib ) and ( 11c ) for 
the cumulants ar e th e same as the classical equations (Uq), excep t the for t he appearance of the nonlinear function 
W(z,z*) in Eq. (lib), which makes the cumulant equations (lib) and (11c) nonlinear. However, this difference may 
be overcomed by introducing the new variables 

1 



B = N{6a*Sa) 



C = N{(SaY), 



where B is real and C is complex. Using (^), the quantum equations of motion 
( |l9| ) in the following form 

1 



N 



Q(z,z*,B,C, C*), 



(19) 

Lll) may be rewritten in variables 



(20a) 



Q(z,z*,B,C, C* 



1 ( d 2 V 

2 \da 2 



C 



d 2 V 
da* 2 



C* 



d 2 V 
da*da 



B 



C = 2 |'£ 
da 



iB 



[dV* 



\ da 



C + 2 



C 



dV 

da* 

dV 

da* 



B. 



(20b) 



(20c) 



(20d) 



and the corresponding equation for C*(t) that could be obtained from (20c) by complex conjugation. Initial conditions 
for the system (20) are 



B(0) = l/2, C(0) = C*(0) = 0. 



(21) 



Now all variables are of the order of unity and the dependence on t he sm all parameter 1 /N is present explicitly only 
in the expression for the quantum correction to classical motion in (20a). 

Compare the set of classical Hamiltonian equations (13) and equations of motion (|l5|) for classical linear fluctuations 



with the equations of motion (20) for mean values and quantum cumulants. It is evident that both sets of equations 
have the same structure with the following principal differences. First, the quantum equations ( |20cj ) and ( p0d| ) 
for the cumulants, in contrast to the clas sical e quat ions (fl5|), are calculated near the mean value z and take into 
account the quantum correction Q [Eqs. (20a), ( 20b )]. Second, it should be noticed that it is impossible to obtain 
the initial conditions (|2l]) for C and B from the initial conditions for the classical equations (O). However, if one 
considers the case of large N, when the quantum correction is small, the quantum equations (|20|) are identical to 
the classical equations ( |7q ) used in the definition of the maximum Lyapunov exponent. The same conclusion on the 
equivalence of the equations of motion for low-order cumulants and equations of motion arising in the definition of 
the Lyapunov exponent has been obtained earlier for the generalized Tavis-Cammings model in 28 and for systems 
with Hamiltonians consisting of the sum of kinetic and potential energies in . 

The self-consistent system of equations ( p0| ) completely describes the dynamics of quantum fluctuations in the first 
order of 1/N. We shall use these equations for the description of the dynamics of light squeezing at the transition to 
quantum chaos in section 4 after a discussion of the validity of our approach and its comparison with other semiclassical 
methods. 
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III. CONDITIONS OF VALIDITY AND COMPARISON WITH OTHER APPROACHES 



Here we discuss the condition of validity of the 1/iV-expansion (||) and our equations of motion ( p0| ) for the mean 
values and cumulants. In analogy to the classical distance (|l7|), we introduce the "distance" D q for the quantum case 
as 

D,(t) = ll?V2( t ) ~ [<|M 2 )] 1/2 - (22) 

It is easy to see that D q coincides with the "convergence radius" r = [Ke(5a)] 2 + [Im(<5a)] 2 of the 1/iV-expansion (|^). 
Initially D q (0) ~ 1/N <C 1, and if D q (t) <g; 1 during some time interval, then the 1/iV-expansion is well-defined in 
this time interval. 

When quantum correction in Eq. ( |20a| ) is small, equations ( |20| ) are identical to the classical equations ( jl5| ) arising 



in the definition of the classical Lyapunov exponent (18). Thus, for classically chaotic motion, we have exponential 
growth of B(t), C(t), and D q (t) as 

D q (t) a D q {0) exp(Aa;ot), (23) 

where A is the Lyapunov exponent and luq is the characteristic frequency in the classical dynamic system. For example, 
luq = g in the model of the Kerr oscillator From the conditions D q {t) <C 1 and D q (0) ~ 1/N, we have the following 
estimate for the time scale of validity of our semiclassical approach 

i<t* = — InJV. (24) 
Xujo 

In contrast, for classically regular motion D q (t) increases with time only power wise: D q (t) ~ t 1 , where the index 
7 ~ 1 depends on the system under study. For example, 7=1 for nonlinear oscillator with the Hamiltonian (||) for 
e = 0, as well as for integrable nonlinear oscillators with power wise higher-order nonlinearity jh|. As a result, the 
time scale of validity of our approach in the case of regular motion is 

t < t* = oj ~ l N th (25) 

and it is much greater than ( p^ ) for chaotic motion. The time scale t* has very simple physical meaning: during time 
interval t* we have localized wave packet with center moving along the classical path. 

The scale (|24|) of the well-defined quantum-classical correspondence for chaotic systems was first introduced in plj , 
and now it is investigated in details in different systems (for a review, see |l6|]). The time scale ( pi] ) is very short 
for small N . However, under the conditions of modern experiments on light squeezing, where the average number of 
photons N involved in the nonlinear interaction is 10 3 — 10 12 , t* is of the order of 10 — 100 of the system's characteristic 
periods luq 1 and thus it looks quite reasonable. 



uantum fluctuations with that widely used in 
8pQ-E2|. This approximation also assumes 



We now compare our approach to the description of the dynamics of qi 
quantum optics and called the generalized Gaussian approximation [p3 ,| 
the existence of a well-localized, almost Gaussian wave packet throughout the time of quantum evolution and in this 
respect it is very closed to our consideration. More precisely, the generalized Gaussian approximation consists in an 
assumption that the Fourier transform of the quantum distribution function, i.e. the quantum characteristic function, 
is Gaussian for any moment of time |p8| . Such a quantum state corresponds to the superposition of the coherent signal 



and small quantum noise |23|. For such a state, only the first- and second-order cumulants are nonzero. It possesses 
the expression for higher-order cumulants in terms of only the first and second-order cumulants and truncation of an 
infinite hierarchy dynamic system for cumulants is possible. 

We compare our system ( p0[ ) with the dynamic equations for the mean values and cumulants obtained within 
the generalized Gaussian approximation for several popular models of quantum optics: second harmonic generation 
described by the Hamiltonian 

H = ujib\bi + uj 2 b\b2 - {rd)\b\ + H.c.^ , [6,6+] = 1, uj 2 = 2u u (26) 

the problem of nondegenerate optical 3-wave mixing described by the Hamiltonian 

3 

H = ^ujjb]^ ■- (k&i^ + H.c.) , w 3 =wi+w 2 , (27) 



G 



and for the periodically forced nonlinear oscillator with the Hamiltonian (0) . Considerations of these problems 
within the generalized Gaussian approximation are presented in Refs. [40|, flUf , and jl2|], respectively. For the forced 
nonlinear oscillator, we found that our self-consistent set of equations (11) and ( |l2| ) coincides with the corresponding 
basic equations of E2] up to terms of the order of 1/N 2 , and for the problems of nondegenerate and degenerate optical 
3-wave mixing, our approach gives equations that are identical to the corresponding equations of motion obtained 
within the generalized Gaussian approximation poj^lp^ ]. 

We turn now to the discussion of squeezing at the transition from regular to chaotic motion. 



IV. SQUEEZING AT THE TRANSITION TO QUANTUM CHAOS 

Define the general field quadrature as X(8) — aexp(-id) + <r exp(i8), where 8 is a local oscillator phase. A state 
is said to be squeezed if there is some phase 8 for which the variance of Xg 

((SX) 2 ) = ((5a) 2 ) exp(-i28) + ((6a*) 2 ) exp(i20) + 2(\Sa\ 2 ) + 1/N (28) 



is less than the variance for coherent state or the vacuum |22|] . In variables (fL9|) , the condition of squeezing takes the 
form 

S(8) = 2B + Ccxp(-i26») + C* exp(i26>) < 1. (29) 

The minimum of the variance of the general field quadrature fl29| ) with respect to its dependence on the local oscillator 
phase 8 is reached at 8 m i n , determined as 

9 mi n = (<p- ?r)/2, (30) 

where ip is the argument of the cumulant C [C = \C\ expfiip)]. For 8 — 8 m i n , the condition of squeezing (|2^) is 

S=S(8 min ) = 2(B-\C\) <1. (31) 

The value S determines the minimum half-axis of the quantum noise ellipse p3[ . The condition (^) is called principal 
squeezing because it gives the maximum squeezing measurable by homodyne detection [[l8| , p3| . 

Let us now compare the dynamics of the principal squeezing for classically regular and chaotic motion. Initially 
the Gaussian wave packet spreads when it propagates through a nonlinear medium. However, there still exists the 
time interval of the well-defined quantum-classical correspondence ( p4[ ) or (p5[) dur ing which the wave packet's center 
follows a path in phase space governed by the semiclassical equation of motion (|20a|). Moreover, because our equations 



of motion for cumulants (20c) and (20d) in fact coincide with the equations arising from the definition of the maximum 
Lyapunov exponent, we can apply simple physical arguments to the strong deformation of the classical phase volume 
at chaos for the prediction of the strong squeezing of the noise ellipse at quantum chaos in the semiclassical limit. 

Due to the presence of the strong (exponential) local instability inherent in the underlying classical chaotic dynamics, 
a quantum noise ellipse may be strongly stretched in one direction and squeezed in another direction. As a result, the 
value of principal squeezing S in Eq. (pi]), which correlates with the minimum half-axis of the quantum noise ellipse, 
on average exponentially decreases in time 

S(t) ~ exp(-\cu t). (32) 

The stretching and squeezing of a noise ellipse at quantum chaos is much stronger than for the case of regular and 
stable dynamics, when the distance between two initially closed trajectories in phase space increases in time in a 
power wise way resulting in only a power wise decrease of the principal squeezing in time 

S(t) ~ (u; t)-P, (33) 

where the constant f3 ~ 1. As an example of such time-dependence, we consider the Kerr oscillator (JsJ) without an 
external field (e = 0). In this case the model is integrable, and from an exact solution in the limit N 1 j49|, we get 

S(t) = 1 + 2\a Q \ 2 gt [\a \ 2 gt - (l + \a^g 2 t 2 f /2 ] 9 ^ (6|a |V) _1 , (34) 

where ao is some initial condition ([ckq | — 1). The same result could be obtained directly from our semiclassical 
approach neglecting quantum correction in Eq. j20| ) and combining it with Eq. ( |3l] ) 46 1. 
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Making a comparison of squeezing for regular and chaotic motion, several comments are necessary. 

1. Both formulas ( |32| ) and ( fffi ) are obtained within the pure classical picture when quantum corrections are neglected. 
The deviations from dependencies (|3|) and are expected when quantum corrections become sufficient, especially 
at t ~ t*. 

2. Squeezing for chaotic dynamics is exponential only on average. Actually, the rate of phase volume deformation 
and, consequently, the rate of quantum noise ellipse deformation and squeezing are directly related in the semiclassical 
limit (N 3> 1) to the degree of local Lyapunov instability D q [Eg. (pffi) ] in the system, which at some moments in 
time may be different from the maximum Lyapunov exponent ( |18| ) measured at asymptotics t — > oo. Such behavior 
is typical, for instance, for a chaotic trajectory that spends some time near a stability island. In this case, the local 
instability and statistical properties of chaotic motion are weak |30| . Finally, the trajectory escapes to a large chaotic 
sea, and the local instability becomes strong. Of course, the Lyapunov exponent is positive in both cases. Formally, 
the resulting nonmonotonic dependence of squeezing on time could be modelled by a slowly varying dependence of 
the parameter A [Eq. (|32|)]. We will say more on this problem in the next section. 

3. Along with wide class of integrable systems with stable regula r dy namics, there is a small class of systems with 
regular but unstable dynamics, for which the distance D q given in (|22j) increases in time exponentially. In such a case 
the systems are in the state near the bifurcation point between different dynamic behaviors and they are exponentially 
unstable. The mechanism of enhanced squeezing in unstable systems with regular dynamics is very similar to that 
discussed at the transition to quantum chaos and has been discussed in detail in p^-^7t . 

4. Considering squeezing at chaos, we have not discussed yet the influence of one of the main characteristics of chaos: 
folding of the phase volume that in addition to stretching, is present in any bounded Hamiltonian system. For strong 
chaos, strong and multiple folding of the phase volume appears during the time scale of well-defined quantum-classical 
correspondence. The perimeter of the phase volume increases in time exponentially, and eventually the phase volume 
envelope gets some fractal-like structure Moreover, the final shape of the phase volume is very sensitive to 
small changes of the initial conditions and/or parameters. The complex and unstable evolution of the phase volume 
is displayed in the time-dependence of squeezing |28|. For strong chaos and long enough time, squeezing becomes 
strongly dependent on tiny variations of the system's parameters. Such a regime of squeezing was called in psll 
unstable squeezing. For unstable squeezing, the range of the local oscillator phase, for which squeezing is possible, 
becomes so narrow that it makes observation of squeezing practically impossible j2S| . 

It should be noticed that at the same time when multiple folding of wave packet appears, the quantum-classical 
correspondence breaks down. For strong chaos this time scale is very short. In contrast, for weak chaos, the time 
scale of the well-defined quantum-classical correspondence and time scale of strong folding may be rather long in 
the semiclassical limit resulting in stable and enchanted squeezing. The same arguments should be valid for regular 
trajectories located near the chaotic motion in phase space. Actually, during a finite time interval stability properties 
of a given trajectory are determined by the time-dependence of distances ( |l7| ) or (^2|) between two initially nearby 
trajectories. At a finite time the value of D q for the regular trajectory is often practically indistinguishable from the 
nearby chaotic trajectory |51J resulting in the same degree of squeezing. 

We will illustrate a general picture of squeezing at the transition to quantum chaos through an example of the 
model of the kicked nonlinear oscillator in the next section. 

Concluding this section, let us compare qualitatively the maximum achievable degree of principal squeezing for 
regular and chaotic motion during the time interval of the well-defined quantum-classical correspondence. First 
consider the case of squeezing at chaos. Principal squeezing on average is exponential (|32|), and at the end of time 
interval ( |24| ) it should be not less than S(t*) ~ exp(— \u>ot*) — l/N. For regular dynamics, the principal squeezing 
has a power wise time-dependence ( |33"|) and it should be not less than S(t*) ~ (woi*)~^ — N~^^ , where constants 
7 ~ 1 and /3 ~ 1, and we used the estimate ( p5| ) of t* for regular motion. 

Thus these rough estimates show that the degree of squeezing at time t* is comparable in the cases of regular and 
chaotic dynamics, but because t* is much shorter for chaotic dynamics than for regular, squeezing is much faster in 
the chaotic systems. 



V. EXAMPLE: KICKED NONLINEAR OSCILLATOR 

Consider a nonlinear oscillator interacting with a time-periodic field as given in Eqs. (|) and (|). This model 
describes, for example, a high-Q cavity filled by a medium with Kerr nonlinearity and excited by an external laser 
field . The same effective Hamiltonian may also govern the interaction of a laser field with a high-density excitons 
in a semiconductors f32l. The different variations of model discussed are very popular in both quantum optics ]53| ] 
and quantum chaos |nj42|,Q studies. 

In this paper we choose the form of F(t) in (|5|) as a periodic sequence of kicks 
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S T (t)= J2 S(t-nT), 



(35) 



where 5(t) is the Dirac (5-function. In an experiment, a sequence of short light pulses can be generated by a mode 
locked laser. The use of the sequence of kicks possesses to obtain discrete maps instead of differential equations and 
it sufficiently simplifies computations and reduces numerical errors, which is especially important when we analyze 
the influence of small quantum corrections (N 3> 1) on the dynamics of squeezing. 

In what follows we shall use scaled variables t' = gt and A' = A/g measuring time and detuning in the units of the 
coupling constant g, and then we omit primes, which is formally equivalent to choosing g = 1 in the Hamiltonian (|^). 

The classical equation of motion ([l^) obtained from the Hamiltonian (||) has the form 



dz 

i— = Az + g\z\ 2 z + sS T {t). 

From Eq. (|36] ) and using a standard technique (l^j3l[], we have the classical map 

z n+1 = F(z n ), F{z n ) = exp[-£T(A + \z n - ie\ 2 )]{z n - e), 



(36) 



(37) 



where z n is the value of z just before n-th kick, z n = z(nT — 0). It may be shown analytically |54|| that the approximate 
criterion of transition from regular motion to strong and global chaos for the map (B7|) in the limiting case e <C 1 is 



K = 2eT > 1. 



(38) 



We illustrate the dynamics governed by the map (|3^) in Fig. |l|, where phase portraits (a,c,e) and the time dependence 
of intensity \z n \ 2 (b,d,f) are shown for the cases of regular motion, weak chaos, and strong chaos correspondingly. 

We now turn to the dynamics of quantum cumulants and the mean values with quantum corrections. Using the 
smallness of 1/N, we obtain from the semiclassical equations of motion ( pp| ) and the classical map d37| ) the following 
coupled maps describing the dynamics of mean values and cumulants (details of the derivation of these maps are 
presented in the Appendix) 



z n+ i = F(z n ) + jyQn, 



(39a) 



C n+ i = exp [-i2T (A + \z n - ie\ 2 )] x 

j-2 {z n - ief (T 2 \z n - ie\ 2 + iT) B n + (l - 2iT\z n - ie\ 2 ) C n 

- T 2 (z n - ief \{z n - ief C* + c.c. 



(39b) 



B n+1 = (2T 2 \z n - ie\ 4 + 1) B n + T 2 \z n - ie\ 2 {z n - ief C* + c.c 



— iT (z n — ie) C* + c.c. 



(39c) 



where F(z n ) in (39a) is the classical map introduced in (|37j) and the quantum correction Q n in ( p9a| ) has the form 

Q n = A 1 (z n ,z*)C n +A 2 (z n ,z*)C*+A 3 {z n ,z^)(B n - 1/2), (40a) 



A 1 {z n , z* n ) = _ e -^( A +l^— I 2 ) [-iT(z* + fe)(2 - %T\z n - ie\ 2 )] , 



A 2 (z n , zl) = L e -mA + \^) ^ T 2 [zn _ fe) 3] ( 



A 3 (z n , z*) = -e-^ A+ ^-^T(z n - te)(T\z n - ie\ 2 + 2i). 



(40b) 
(40c) 
(40d) 
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Following (21) the initial conditions for the maps ( p9[ ) are Bq = 1/2, Cq = 0, and arbitrary zq of order of unity. 



The self-consistent set of maps (39) and (|4Cj) determines the dynamics of the quantum fluctuations for kicked 
nonlinear oscillator in first order of X/N. 

Now we want to compare the time-dependence of the principal squeezing S in (|3l|), the degree of local instability 
(p2|), and the quantum correction Q in ( fl0| ) for three characteristic cases of classical dynamics: regular motion, mild 
chaos and hard chaos. To avoid the dependence of D q on N it is useful to introduce new normalized distance 

d = ND q . 

Phase portraits and time-dependencies of the intensity \z\ 2 for the three characteristic cases of classical dynamics 
governed by the map ( |37| ) are shown in Fig. 1. The time-dependences of log 10 S', log 10 d, and log 10 Q for a large but 
finite number of quanta N — 10 9 are shown in Fig. 2, where curves 1, 2, and 3 correspond to the cases of regular 
motion (A = 0) and chaotic motion with an increasing value of the Lyapunov exponent A, respectively. As is evident 
from comparison of Fig. 2a and Fig. 2b, the strongest local instability determines the highest degree of squeezing. It 
should be mentioned that the difference in the magnitude of principal squeezing for chaotic motion and regular motion 
is of several orders during only several kicks. A power wise time-dependence of d (Fig. 2b, curve 1) is assisted by 
the corresponding slow growth in time of the quantum correction (Fig. 2c, curve 1). In contrast, for chaotic motion, 
the growth of d and Q is exponential (curves 2 and 3 in Fig. 2b and 2c) resulting in a logarithmic dependence of the 
applicability of semiclassical approach on the number of quanta in J24|). 

The increment of local instability and the corresponding degree of squeezing are complex functions of the initial 
conditions or parameters. As a result, the dependence of squeezing on initial conditions or parameters may be rather 
complex. In Fig. 3, we plot the minimum value minlog 10 5' of the principal squeezing during seven kicks and degree 
of the local instability log 10 G? as a function of the kick's period T at fixed initial conditions, e — 0.1, and number 
of quanta N = 10 9 . The regions of parameter values with mainly regular or chaotic motion in the classical limit 
(N — > oo ) are marked in the center of the plot. The first large region of chaos begins at T cr ks 5.4, which is in good 
agreement with the approximate criterion ( ^8|) of the transition to chaos. Of course, there are some singular values 
of the parameter T corresponding to regular motion in the range of T that is mainly filled by chaotic trajectories. 

As is evident from Fig. 3, there is a direct correlation between the degree of squeezing and degree of instability 
log 10 rf. Because we present the dependences of log 10 d and squeezing on T for only seven kicks, the degree of instability 
and squeezing may be comparable for some points from regular and chaotic ranges of T. However, generally, as it 



follows from (38), the degree of instability increases with increasing T and correspondingly we have an increase of 
squeezing (Fig. 3). 

In Fig. 4 we plot log 10 cZ and the minimum of log 10 <S during nine kicks as a function of another parameter e and 
at fixed T and initial conditions. Again, the correlation of local instability and squeezing as well as enhancing of 
squeezing at the transition to chaos are visible. The same dependences but for three kicks instead of nine are shown 
in Fig. 5a. In spite of the correlation between the degree of squeezing and log 10 e? it is also evident in Fig. 5a that the 
increase of squeezing at the transition from regular to chaotic motion (e cr « 0.5) is smoothed more than in the case 
of nine kicks shown in Fig. 4. 

As we already mentioned, the degree of local instability may itself depend on time. This is illustrated in Fig. 5b, 
where in an interval of 3 — 8 kicks the instability for e = 1.8 (curve 1) is greater than for e = 2.22 (curve 2). As a 
result, the squeezing during 3 kicks at e = 1.8 (point 1 in Fig. 5a) is greater than for e — 2.22 (point 2 in Fig. 5a). 
After eight kicks we have inverse picture: the instability for e = 2.22 is greater than at e = 1.8 (see Fig. 5b) resulting 
in enhanced squeezing for the parameter's values of point 1 in comparison to squeezing for the parameter's values of 
point 2. 

Summarizing our findings from Figs. 3-5, we see that in the semiclassical limit the degree of squeezing has a one-to- 
one correspondence to the degree of local instability in the system. For a short time of the order of several kicks, the 
dependence of squeezing on the Lyapunov exponent, which is calculated in the limit t — > oo, is not well pronounced. 
Moreover, squeezing for regular motion located near the border of the transition to chaos is comparable to squeezing 
for weak chaos. The dependence of squeezing on the value of the Lyapunov exponent becomes more pronounced for 
moderate times of the order of ten kicks. 

We compare now the time-dependence of principal squeezing for different N. The dependence of S on t at N — > oo 
(curve 1) and at N = 10 7 (curve 2) is presented in Fig. 6 for the case of mild chaos. Curves 1 and 2 coincide during 
seven kicks, which determines the time interval of well-defined quantum-classical correspondence for such a choice of 
parameters. Moreover, our numerical calculations demonstrate that if the number of quanta increases in two orders 
up to A = 10 9 , the squeezing is indistinguishable from squeezing at N — > oo up to fifteen kicks. This finding is well 
explained by remembering In N dependence ( pi]) of the time scale. 

Thus we demonstrate that an increase in the average number of photons N results in corresponding increase of 
the time interval for the applicability of our description of squeezing dynamics. On the other hand, the number of 
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pho tons N > 10 7 initially pumped to the system is quite realistic for contemporary experiments on light squeezing 

MM- 

Let us now discuss the stability of chaotic squeezing. The time-dependence of the optimal local oscillator phase 
Omin given in ([30j) for the cases of regular and mild chaos is shown in Fig. 7 for two slightly different initial conditions 
Zo = 1 and zq = 1.001. The deviation of the optimal local oscillator phase values related to different zq in the case 
of regular dynamics is small during all time evolution presented in Fig. 7 (see dashed and dotted curves). For chaos 
(T = 7, e = 0.1), the deviation is sufficient only at the fifth kick (boxes and crosses in Fig. 7). Thus, for not very 
strong chaos, there is a time interval during which squeezing is enhanced and still stable. 



VI. CONCLUSION 

In Ref. a fundamental question was asked: "What is the best squeezing device?" The answer founded in |2|] 
is based on the fact that "the quantum noise reduction is effective below and above the bifurcation points between 
dynamical regimes" . The reason for enhanced squeezing consists in the fact that "the same classical linearized 
equations are used to study the stability of the system and to calculate its quantum fluctuations " p5j . It should be 
stressed that only integrable or near integrable systems with regular dynamics were discussed in |25|. The results of 
our paper demonstrate that in addition to a very narrow class of integrable systems near the threshold of instability, 
there is another wide class of potentially "the best squeezing devices" : systems at the transition to quantum chaos 
operating during the time interval of the well-defined quantum-classical correspondence. 

In this paper we have presented a semiclassical theory for the quantum systems with one and half degrees of freedom 
and illustrated our findings on the model of kicked nonlinear oscillator. Direct numerical observation of enhanced 
squeezing in the model of a kicked quantum rotator at the transition to quantum chaos (nonperturbative approach) 
will be presented elsewhere fl5a]. 
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APPENDIX A: 



In this appendix we briefly describe the derivation of the coupled maps (|39|), and (|4(j) for the mean values and 
second order cumulants. 

Begin with the equations for mean values ( |39a| ) . Substituting the Hamiltonian of the nonlinear oscillator (||) with 
F(t) given in ( |35| ) to (20a) and using the smallness of quantum correction N~ 1 Q, we first neglect the influence 
of the quantum correction and find the solution of the classical equations of motion given by the map ( |37| ) and 
then substitute this classical solution into the expression for the quantum correction Q(z — ► z c i,z* — > z* t ). Such a 
perturbative approach is valid during the time interval ( p4| ) of the well-defined quantum-classical correspondence for 
chaotic systems, and the time interval ( p5[ ) for the systems with regular dynamics. 

To find the explicit form of the quantum correction Q, we remember that it has the form of the second differential 
of z(t) (see ( |l2|) and (|20b|) ). Determining d?z n +i , we obtai n the expression ( |40| ) for the quantum correction Q 



We now turn to the derivation of the maps (39b) and (39c) for cumulants B and C. As it has been shown in 
Sec. 2, the equations of motion for B and C may be obtained by linearization of classical motion equations near the 
classical trajectory taking into account the quantum correction, and by using substitution B — > \dz\ 2 , C — > (dz) 2 . We 
adopt this finding for the derivation of maps for B and C . Utilizing again the perturbative approach, we first find 
the differential of the map (39a) neglecting the influence of the quantum correction 



DZr, 



e - lT ( A+|z "- l£|2 ){-iT(z„ - is) [{z n - is)dz* + « + ie)dz n ] + dz n }. 



(Al) 



From (|A1|) we obtain the maps for quadratic variables B n = \dz.g \ 2 an d C n = {dz n ) 2 in the forms (39b) and (39c) 



Finally, we suppose that the variables z n , z* included in Eqs. (39b) and (39c) are calculated with account of quantum 
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correction, i. e., using formula ( p9a| ). The perturbative approach used in the derivation of ( p9q ) and (39c) is valid 
during time interval of well-defined quantum-classical correspondence when the quantum corrections to the classical 
equations are small. To avoid misprints we checked the expressions for our coupled maps ( |39| ) and ( J40| ) by symbolic 
computation in the package MATHEMATICA. 
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FIG. 1. Nonlinear dynamics of kicked classical nonlinear oscillator (N — > co): phase portrait and time dependence of the 
intensity \z n \ 2 for the regular behavior at T = 3 (a, b), for weak chaos at T — 6.4 (c, d), and for hard chaos at T = 10 (e, f). 
The initial condition is zq = 1, and e = 0.1, A = 1. Time is measured in the number of kicks. 



FIG. 2. Time dependence of the principal squeezing (a), of the distance between two initially closed trajectories (b), and of 
the quantum correction to the motion of wave packet center (c). Curve 1 corresponds to the regular dynamics at T = 3, curve 
2 to the mild chaos (T = 7), and curve 3 to the hard chaos (T = 10). The average photon number is N = 10 9 , The initial 
condition and other parameters are the same as in Fig. 1. 



FIG. 3. In the upper part of the figure, the logarithm of the distance between two initially closed classical trajectories 
log 10 d after seven kicks is plotted as a function of kicking period T defined in dimensionless units (see main text). In the lower 
part, the minimum of the logarithm of the principal squeezing minlog 10 5' during seven kicks as a function of the kick period 
T is plotted for the average photon number N = 10 9 . In the center of the figure, big intervals of T with primary regular and 
primary chaotic behavior are marked. The amplitude of perturbation and detuning are fixed: e = 0.1, A = 1. The initial 
condition is the same as in Fig.l. 
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FIG. 4. Local Lyapunov rate log 10 d after nine kicks and the minimum of principal squeezing during nine kicks as functions 
of dimensionless perturbation amplitude e at fixed T = 1. Other parameters are the same as in Fig. 3. 

FIG. 5. (a) Same as in Fig. 4 but after only three kicks instead of nine. Point 1 corresponds to e = 1.8 and point 2 to 
e — 2.22. (b) Time dependence of the classical local Lyapunov rate. Curves 1 and 2 are plotted for the same values of e as in 
(a) and at T = 1. 

FIG. 6. Time dependence of the principal squeezing for chaotic motion (T = 7, e — 0.1). Curve 1 corresponds to the classical 
limit iV — ► co, and curve 2 to N = 10 7 . 

FIG. 7. Time dependence of the optimal local oscillator phase min in the case of mild chaos at T = 7, e = 0.1, and for 
different initial conditions: zo = 1 (boxes), z = 1.001 (crosses). Compare with time dependence of 6 min for regular motion at 
T — 3 and e = 0.1 showed by dashed line for zo = 1 and by dotted line for zo = 1.001. 
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